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Universe was a blackbody spectrum in vacuum. We derive the electromagnetic 



spectrum, based on the Fluctuation-Dissipation Theorem that describes the 



electromagnetic fluctuations in a plasma. Our description includes thermal 



differs from the blackbody spectrum in vacuum at low frequencies. In partic- 



and collisional effects in a plasma. The electromagnetic spectrum obtained 
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ular, concentrating on the primordial nucleosynthesis era, it has more energy 
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for frequencies less than 3 to GcOpe, where ujpe is the electron plasma frequency. 
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I. INTRODUCTION 



It is usually assumed in Cosmology that the primordial plasma was a homogeneous 
plasma and that the electromagnetic field was a blackbody spectrum in vacuum. These 
assumptions are used, for example, in standard Big Bang Nucleosynthesis calculations. 
Deviations from a blackbody spectrum in vacuum and a homogeneous plasma can affect 
Primordial Nucleosynthesis and in this letter we concentrate on this epoch. 

In the epoch of Big Bang Nucleosynthesis, the universe was a thermal bath of photons, 
electrons, positrons, baryons and neutrinos. The usual treatment considers the universe as 
a homogeneous plasma in thermal equilibrium and the electromagnetic field as a blackbody 
spectrum in vacuum. For example, in the energy density calculation, the energy density of 
the photons is given by the energy density of the blackbody spectrum in vacuum. 

In this manner, the Primordial Universe is treated as an ideal gas: collective effects are 
assumed to be negligible. However, a plasma differs from an ideal gas due to correlations, 
for example, two-particle correlations. In an ideal gas, the distribution function for two- 
particles is: -F2(xi, X2) = Fi{xi)F2{x2). In a plasma we have an additional term F2{xi, X2) = 
[1 + Pi2{xi,X2)]Fi{xi)F2{x2), where Pu is the Debye-Huckel screening. (P12 oc exp{— \ 
X2 — Xi I /Xd)/ \ ^2 — Xi I, where is the Debye length.) 

A plasma, even in thermal equilibrium, has fluctuations, that is, the physical vari- 
ables such as temperature, density and electromagnetic flelds, fluctuate. Even for a non- 
magnetized plasma, where the average magnetic fleld is zero, (B) = 0, the squared average 
is not zero: (5^) 7^ 0. 

The study of the electromagnetic fluctuations in a plasma has been made in numerous 
studies, including those of Dawson, Rostoker et al., Sitenko and Gurin and Akhiezer et al. 
Most of the results are compiled in Sitenko and Akhiezer et al. 0. 

The electromagnetic fluctuations are described by the Fluctuation- Dissipation Theorem 
0. The intensity of such fluctuations is highly dependent on how the plasma is described, 
for example, on the dissipation mechanisms present in it. It is necessary to describe the 



2 



plasma in the most complete way. 

Cable and Tajima (see also 0,0) studied the magnetic field fiuctuations, for several 
cases. Two of there descriptions concern the primordial plasma, which is an isotropic, non- 
magnetized and non-degenerate plasma. In particular, they studied: a) a cold, gaseous 
plasma and b) a warm, gaseous plasma described by kinetic theory. 

In their study. Cable and Tajima [^j in case (a) used the cold plasma description with a 
constant collision frequency. In case (b) they analyzed the spectrum only for low frequencies, 
with the warm plasma description for phase velocity uo/k less or equal to the thermal velocity 
of the electrons, fe, and the ions, fj, in a collisionless description. 

Through the study of the electromagnetic fluctuations, that is, the magnetic and electric 
field transverse fluctuations, given for example by the Fluctuation-Dissipation Theorem, a 
reliable electromagnetic spectrum can be obtained. In this letter, we present a model that 
extends the work of Cable and Tajima 0. It includes in the same description collisional 
and thermal effects. By using the Fluctuation-Dissipation Theorem relations, we derive the 
electromagnetic spectrum in a plasma. We concentrated, in particular, in the primordial 
nucleosynthesis era, in the case of the electron-positron plasma at high temperatures and the 
electron-proton plasma at low temperatures. The first case is the plasma at the beginning of 
the primordial nucleosynthesis, when the number of electrons and positrons was comparable 
to the number of photons. The second case is the plasma at low temperatures, after the 
annihilation of the electrons and positrons where the number of electrons and protons was 
Tj ~ 10^^*^ smaller than the number of photons. 

In Section II we present the expression of the electromagnetic fiuctuations, in Section III 
our model, and in Section IV a discussion of our results and conclusions. 

II. ELECTROMAGNETIC FLUCTUATIONS 

The spectrae of fiuctuations of the magnetic and electric transverse fields, for an isotropic 
plasma, given by the Fluctuation- Dissipation Theorem (for the deduction see Sitenko and 
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Akhiezer et al. are 



where Et is the transverse dielectric permittivity of the plasma. 

An intuitive way to understand the above expressions, is that the Fluctuation-Dissipation 
Theorem takes into account the emission and absorption processes in a plasma, and knowing 
that in equilibrium they are equal, the fluctuation level is obtained. 

To determine the fluctuations of the magnetic and electric transverse fields in a plasma in 
equilibrium or quasi-equilibrium, it is sufficient to know the transverse dielectric permittivity 
of the plasma, in particular, the dissipation mechanisms present for each frequency and wave 
number (Jm Et)- This depends on the treatment used to describe the plasma. 

Another important feature of the electromagnetic spectral distributions can be seen from 
Eq. (1). The equation Et{'jJ-,'^) — = determines the transverse eigenfrequencies of 

the plasma. In the transparency region (Jm et <^ Re £t), the magnetic and the electric 
field transverse spectrae have 5-function-like maximae near the eigenfrequencies (i.e., the 
frequency spectrum of the fluctuations contains only the transverse eigenfrequencies in the 
plasma). Knowing that photons have the dispersion relation in the plasma u"^ = + k'^c^, 
we note that for frequencies uj ujp the eigenfrequencies dominate and the electromagnetic 
transverse field spectrum behaves like a blackbody spectrum. 



III. BGK COLLISION TERM 

To describe completely the plasma, we present a model that includes thermal effects 
as well as collisional effects. For this we need a kinetic description that takes into account 
collisions. We use the Vlasov equation in first order |]^ (in the plasma parameter g = 
where A^) is the Debye length and n is the particle density). The collisions are described by 
the collision term in this equation. 

Charged particles simultaneously interact with all the particles in the Debye sphere, 
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which is a large number and the Fokker-Planck colhsion term is most appropriate for a fully 
ionized plasma. It describes the effect of the microscopic fields produced by all the particles 
in the plasma. However, the kinetic equations with the Fokker-Planck collision term are very 
hard to solve. Therefore, we used the BGK collision term as a rough guide to the inclusion 
of collisions in the plasma. The BGK collision term describes the binary collisions as used 
by Cable and Tajima |^ and is given by, (for a derivation see Clemnow and Dougherty or 
Alexandrov et al. p) 



dt 



= -Vif - fmax) , (2) 

C 



where rj is the collision frequency (considered constant and equal to the Coulomb collision 
frequency), fmax is given by fmax{^,'^,t) = A^(x, t)/o(v)/A^o, where the number density is 
N = Nq + Ni and / = /o + /i, /o being the unperturbed Maxwellian distribution. (This 
collision term conserves the number of particles.) Substituting this collision term in the 
equation of Vlasov in first order, it is straightforward to obtain, for an isotropic plasma, the 
transverse dielectric permittivity Et, (generalized for several species): 

M-,k) = l + E^f^Vf^l , (3) 

where a is the label for each species of the plasma, Va the thermal velocity for each species 
and Z{z) the Fried & Conte function 0, 

1 (4) 

\/7r J-oo t — z 



(If relativistic temperature effects are included, the substitution Upa — > cUpa/i/r is made.) 

It is interesting to comment and emphasize that both the cold plasma description and the 
warm plasma coUisionless description, used by Cable and Tajima Q], are particular solutions 
of this model. For | z |^ ^ 1, where z = {uj + irj)/ ^/2fcfe, we obtain the cold plasma dielectric 
permittivity and for t] —>■ the warm plasma coUisionless dielectric permittivity. 

We substitute the transverse dielectric permittivity in Eq. (1), obtaining the magnetic 
and the electric transverse field spectrum, (5^)^^^ and {Et^)^^^ respectively. Integrating 



over wave number (and dividing by (27r)^) gives us the magnetic and electric field transverse 
spectrae, {B'^)^ and {Enp'^)^, respectively. 

Our model uses a kinetic theory description with a collision term that describes the 
binary collisions in the plasma and a cut-off has to be taken, since for very small distances 
the energy of the Coulomb interactions of the particles exceeds their kinetic energy which 
violates the applicability of the condition of the perturbation expansion (in the plasma 
parameter g <^ 1). This occurs approximately for distances r^m ~ e^/T, or more exactly, 
the distance of closest approach between a test particle and an electron in a plasma, r^m = 
{kmaxy ^ — [Mmv'^/ (m + M) | eg |] , where M, v and q are respectively, the mass, velocity 
and charge of the test particle 0. 

Treating properly the effects of distant encounters, the cut-off procedure can be removed. 
This was shown by Thompson and Hubbard, and Hubbard, in several works [9], analyzing 
the Fokker-Planck equation and its coefficients. They showed that the cut-off procedure is 
unnecessary, when higher order terms in the Fokker-Planck equations are maintained. The 
kinetic equations with the Fokker-Planck collision term are very hard to solve, however. We 
used the BGK collision term where a cut-off is necessary. We took k^ax consistent with this 
collision term. 

The electromagnetic spectrum is obtained by summing the magnetic and electric field 
transverse spectrae: 

S{u;)=^-^+^-^. (5) 

In Figure la, we plot the electromagnetic spectrum S{uj) given by Eq. |^ (divided by 
the normahzation So = UpJisT / c") vs u/upe- We study an electron-positron plasma at 
T = 10^° K and rig = 1.4 x 10^^ cm'^. This is the plasma at the beginning of the Pri- 
mordial Nucleosynthesis era, when neutrinos decoupled from the plasma and the neutron 
to proton ratio become frozen-in. (This ratio essentially determines the primordial "^He 
abundance.) The dotted curve is S{uj) (our model) compared to the blackbody spectrum 
in vacuum (the solid curve). In Figure lb, we plot the same curves as in Figure la, but 
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extended to high frequencies. In Figure Ic and Id, we did the same as Figure la and lb, 
but for an electron-proton plasma with T = 10^ K and = 5.4 x 10^^. In the epoch of 
Primordial Nucleosynthesis, at lower temperatures, the electrons and positrons annihilate 
and the plasma is reduced to a plasma of protons and electrons. 

It can be seen that we obtain the blackbody spectrum naturally for high frequencies, 
and in the case of the high-temperature plasma (T = 10^° K), for frequencies a; < 3 ujpe-i the 
spectrum has more energy than the blackbody spectrum. In the case of the low- temperature 
plasma (T = 10^ K) it has more energy for frequencies u < ^Upe- 

IV. DISCUSSION AND CONCLUSIONS 

The unique manner to obtain the electromagnetic transverse spectrum is analyzing the 
magnetic and electric field transverse fluctuations. This is the only manner to obtain in- 
formation, not only about modes that propagate, like photons, but also modes that do not 
propagate. These modes appear, not only at low frequencies but also at high frequencies, 
resulting from the correlations in the plasma. Only at very high frequencies, the photons 
contribute uniquely to the magnetic and electric field transverse spectrae. 

We present a model that incorporates, in the same description, the thermal and coUisional 
effects and used the Fluctuation-Dissipation Theorem that describes the electromagnetic 
fluctuations. We use the Vlasov equation with the BGK collision term. 

The final electromagnetic spectrum for the primordial plasma at the epoch of Big Bang 
Nucleosynthesis behaves like a blackbody spectrum in vacuum for high frequencies. How- 
ever, for low frequencies, it is distorted. It has more energy than the blackbody spectrum 
in vacuum. In the case of the high-temperature plasma (T = 10^° K), this range is for 
frequencies uj < ScUpe and for the low-temperature plasma (T = 10^ K) the range is for 
LU < 6u!pe, where cUpe is the electron plasma frequency. 

This additional energy is due to the collective modes of the plasma. The reason why 
the collective modes of the plasma can have more energy for cu < cUpg than the photons in 
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vacuum, can be understood as follows. Photons are massless bosons with the dispersion 
relation lJ^ = k'^c^. For the energy interval < uj < cjpe, the wave number interval is 
k = to k = ujpe/ c. A relatively small amount of phase space is involved. For the collective 
motions of the plasma, in general, we have a larger amount of phase space. For example, 
for plasmons with energy uo ~ cUpe, the amount of phase space extends to a maximum k of 
ko = Upe/vT, where vt is the thermal electron velocity, which is greater than ujpe/c for the 
photons. In general, for a given frequency for u < Upe, the greater phase space available 
to the collective modes of the plasma (than that of the photons) implies more energy, or a 
higher spectrum. 

This result (the additional energy that appears in the electromagnetic spectrum com- 
pared to the blackbody energy usually assumed) can affect several fields in Cosmology, in 
particular. Big Bang Nucleosynthesis. This extra energy, that has not been previously taken 
into account, causes the Universe to expand more rapidly at a given temperature. In partic- 
ular, it causes the neutrinos to decouple earlier (at a higher temperature) and the neutron to 
proton ratio to freeze-in at a higher value. (In order to estimate the total additional energy, 
we have to add also, the longitudinal energy due to the longitudinal electric field spectrum. 
For T = 0.8 MeV, we obtain for the additional energy, Ap = l%p^, where = Pbb, the 
energy density of the blackbody photon spectrum in vacuum. A complete study, estimating 
the additional energy for diverse temperatures and densities, is in preparation ||T0[|). 

Another interesting aspect is how the additional energy affects the spectrum of the 
microwave background. The additional energy at zdec (the redshift when the spectrum 
was formed) occurred at frequencies ~ 10"^ t^peafc iyjpeak = 2.8kBT/h). However, non- 
linear effects in plasma can bring the additional energy to higher frequencies. This is very 
interesting and should be investigated in the future. 
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FIGURES 

FIG. 1. The electromagnetic spectrum ln[S{uj)/So] vs oo/upe (where 
S{(jj) = {B'^)^/8tt + {Et'^)^/8it and 5*0 = uOpe^ksT / c" is the normaUzation) for: (a) The elec- 
tron-positron plasma at T = 10^° K and ng = 1.4 x lO^^cm"^ (the dotted curve is our model 
and the solid curve is the blackbody spectrum in vacuum); (b) The same as case (a), extended to 
high frequencies; (c) The same as case (a) but for an electron-proton plasma at T = IQ^K and 
Ue = 5.4 X 10^^; and (d) The same as case (c), extended to high frequencies. 
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